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1.  INTRODUCTION 


1.1  The  Problem  of  Control  in  Decentralized  Organizations 

In  recent  years,  ouch  emphasis  has  been  placed  on  decentralized 
decision-making  in  large  organizations.  This  is  due  to  the  fact  that  very 
sophisticated  decision  aids,  large  amounts  of  decision  information,  and  the 
need  for  quick  response  have  made  it  impossible  to  refer  all  critical 
decisions  to'the  top  levels  of  the  hierarchy.  This  thrust  toward  decen¬ 
tralized  control  has  led  to  several  new  problems,  however.  One  such  issue 
concerns  the  need  for  higher  hierarchical  levels  to  monitor  the  decisions 
made  at  lower  levels,  and,  if  necessary,  intervene  in  the  event  of  unsatis¬ 
factory  performance.  In  general,  three  steps  are  involved  in  any  approach 
to  this  problem:  first,  the  organization  and  its  elements  must  be  mathema¬ 
tically  modeled;  second,  the  associated  information  structure  or  flow  must 
be  determined;  and  final1/,  the  appropriate  monitoring  and  control  strategy 
must  be  selected  for  that  organization.  Much  work  has  been  done  in  the 
first  two  areas,  particularly  in  the  modeling  of  strategic  and  tactical 
military  structures  [1].  But  the  third  area,  at  least  in  the  military 
application,  still  relies  heavily  on  traditional  command  structures, 
channels,  and  leadership  theory  [2]. 

Recently,  several  papers  have  been  published  which  formalize  the 
notions  of  incentives  in  organizations  whose  participants  have  different 
objective  functions,  mainly  due  to  competition  between  members  ([3]  and 
[4]).  An  important  result  of  this  work  is  that  virtual  cooperation 
between  the  members  can  be  induced  by  selecting  an  appropriate  control 
strategy,  while  the  members  continue  to  maximize  their  own  objectives. 
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1.  INTRODUCTION 


1.1  The  Problem  of  Control  in  Decentralized  Organizations 

In  recent  years,  much  emphasis  has  been  placed  on  decentralized 
decision-making  In  large  organizations.  This  Is  due  to  the  fact  that  very 
sophisticated  decision  aids,  large  amounts  of  decision  information,  and  the 
need  for  quick  response  have  made  it  impossible  to  refer  all  critical 
decisions  to  the  top  levels  of  the  hierarchy.  This  thrust  toward  decen¬ 
tralized  control  has  led  to  several  new  problems,  however.  One  such  issue 
concerns  the  need  for  higher  hierarchical  levels  to  monitor  the  decisions 
made  at  lower  levels,  and,  if  necessary,  intervene  in  the  event  of  unsatis¬ 
factory  performance.  In  general,  three  steps  are  involved  in  any  approach 
to  this  problem:  first,  the  organization  and  its  elements  must  be  mathema¬ 
tically  modeled;  second,  the  associated  information  structure  or  flow  must 
be  determined;  and  finally,  the  appropriate  monitoring  and  control  strategy 
must  be  selected  for  that  organization.  Much  work  has  been  done  in  the 
first  two  areas,  particularly  in  the  modeling  of  strategic  and  tactical 
military  structures  [1].  But  the  third  area,  at  least  in  the  military 
application,  still  relies  heavily  on  traditional  command  structures, 
channels,  and  leadership  theory  [2]. 

Recently,  several  papers  have  been  published  which  formalize  the 
notions  of  incentives  in  organizations  whose  participants  have  different 
objective  functions,  mainly  due  to  competition  between  members  ([3]  and 
[4]).  An  important  result  of  this  work  is  that  virtual  cooperation 
between  the  members  can  be  Induced  by  selecting  an  appropriate  control 
strategy,  while  the  members  continue  to  maximize  their  own  objectives. 


1.2  Outline  of  a  Proposed  Solution 


It  is  the  purpose  of  this  thesis  to  demonstrate  how  incentive  theory 
may  be  applied  to  the  problem  of  hierarchical  control  within  the  Tactical 
Air  Control  System  (TACS).  In  Chapter  2,  we  supply  the  necessary  back¬ 
ground  in  incentive  theory,  as  it  applies  to  organizations  and  their 
elements  whose  operations  can  be  modeled  by  objective  functions.  Also  in 
this  Chapter,  we  briefly  define  the  structure  of  a  typical  TACS.  Then,  in 
Chapter  3,  we  provide  a  very  simple  model  of  the  basic  building  block  of 
the  TACS,  the  Modular  Control  Element  (MCE).  Finally,  in  Chapters  4  and  5, 
we  propose  two  methods  for  applying  incentive  control  to  the  leader- 
follower  relationship  between  the  Tactical  Air  Control  Center  (TACC)  and 
the  MCE.  The  first  method,  the  piece-wise  linear  approach,  provides  total 
insensitivity  to  variations  in  an  Internal  parameter  of  the  MCE  model, 
while  the  second  method,  the  quadratic  approach,  provides  only  minimum 
sensitivity  but  is  more  analytically  appealing.  Each  of  these  chapters 
also  contains  a  numerical  example  to  help  illustrate  the  concepts. 


*  v  - 
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2.  BACKGROUND 


2.1  Incentive  Theory  in  Organizations 

The  problem  that  is  addressed  in  incentive  control  is  that  of  ensuring 
that  the  members  of  a  particular  organization  behave  in  a  manner  that  opti¬ 
mizes  the  performance  of  the  entire  organization  [3].  We  will  assume  that 
the  participants  in  the  organization  each  have  an  individual  objective 
function  which  must  be  optimized,  denoted  J^.  There  is  also  an  overall 
organizational  objective  function,  Jq,  which  must  be  optimized  by  the 
organizational  leader. 

Suppose  that  each  organizational  element  has  an  objective  function  of 
the  form 


Ji  * 


where  r  is  an  independent  variable,  ^  is  a  characteristic,  internal, 
constant  parameter  of  the  element,  and  n  is  a  parameter  which  depends  upon 
external  data.  Suppose  further  that  the  overall  organizational  objective 
function  is  a  linear  combination  of  the  elemental  objective  functions: 


Jo  "  alJl  +  a2J2  +  +  Vn 


Finally,  we  assume  that  each  is  a  convex  function  in  r,  so  that  each 
element  is  responsible  for  maximizing  its  individual  objective  by  finding 
the  corresponding  r^,  denoted  r^*,  such  that 


Ji  "  f(Si*n,rt*) 
max 

r^*  ■  { r:  dJ^/dr  -  0} 


It  is  apparent  from  Equation  (2)  that  if  each  element  maximizes  its  own  J 
then  will  also  be  maximized. 

The  problem  arises  when  a  participating  element,  which  we  shall  call 
follower,  has  a  different  perception  of  its  own  objective  function  than 
does  the  organizational  manager,  whom  we  denote  leader.  One  interpretation 
of  this  problem  is  that  the  follower  has  a  different  perception  of  its 

internal  parameter,  8^,  although  the  structure  of  its  function  remains 
unchanged.  Suppose  the  leader  assumes  that  the  follower  has  parameter  8', 
while  the  follower  thinks  he  is  characterized  by  a  slightly  different 
constant,  8  (the  subscript  i  is  dropped  from  now  on).  If  the  maximization 
in  Equation  (3)  depends  upon  8,  then  the  follower  will  arrive  at  a  dif¬ 
ferent  r*  than  the  value  calculated  by  the  leader,  denoted  r*'.  Once  the 
leader  observes  this  discrepancy,  he  must  consider  two  issues: 


(1)  Is  the  discrepancy  large  enough  to  worry  about? 

(2)  Is  the  leader's  perception  of  8  more  appropriate? 


If  the  answer  to  both  questions  is  affirmative,  the  leader  may  want  to 
intervene  using  an  incentive  strategy  [3]. 

An  Incentive  strategy  is  a  function,  announced  by  the  leader  to  the 
follower,  which  Induces  the  follower  to  maximize  the  leader's  objective 
function,  while  still  maximizing  its  own  objective.  In  this  case,  it  would 
induce  the  follower  to  behave  as  if  its  internal  parameter  were  8',  instead 
of  8,  thus  maximizing  at  r*'.  In  order  to  implement  any  type  of  incen¬ 
tive  control,  however,  the  leader  must  have  at  least  partial  control  over 
one  of  the  objective  function  parameters,  and  this  parameter  must  sub¬ 
sequently  affect  the  output.  Referring  to  Equation  (1),  we  assume  the 
leader  can  partially  control  the  external  variable  n. 


Now  the  leader  constructs  n  as  a  function  of  r,  the  independent 
variable.  This  function,  denoted  n,  is  announced  to  the  follower,  who  can 
incorporate  n  into  its  objective  function  before  maximizing.  If  n  is 
constructed  properly,  the  follower's  new  objective  function  will  have  its 
maximum  at  r* ' ,  even  though  the  follower  continues  to  use  0,  rather  than 
0'.  This  effect  is  illustrated  in  Figure  1.  In  Figure  1(a),  the  control  n 
is  a  constant  (incentive  not  used)  so  the  follower,  using  0,  maximizes  at 
r*,  instead  of  the  desired  r*'.  In  Figure  1(b),  n  is  a  function  of  r,  so 
constructed  that  the  maximum  now  occurs  at  r*',  the  leader's  maximum 


Figure  1.  Typical  follower  objective  function 


Note  that  J|J(ax  may  not  be  the  same  as  Jmax,  but  this  is  not  important 
since  the  leader  is  assuming  that  0,  which  the  follower  perceives,  is  not 


accurate  anyway.  The  important  fact  is  that  with  the  control,  the  follower 


maximizes  at  r* ' ,  which  corresponds  to  d  ' . 

An  additional  desired  quality  of  this  control  scheme  is  robustness. 
That  is,  the  incentive  strategy  should  be  designed  so  that  the  follower 
will  maximize  at  r* '  for  a  range  of  0,  since  the  leader  will  not,  in 
general,  know  the  exact  8  which  the  follower  will  assume. 

2.2.  The  Tactical  Air  Control  System  (TACS) 

The  TACS  is  a  hierarchical,  command,  control,  and  communications  (C^) 
network,  whose  purpose  is  the  detection,  identification,  tracking  and 
control  of  aircraft  in  a  defined  geographical  area.  Under  the  distributed, 
decision-making  concept,  each  level  of  the  TACS  is  tasked  with  certain  func 
tions,  such  as  surveillance,  and  has  the  distributed  data  and  hardware 
available  to  make  accurate  and  timely  decisions.  The  general  structure  of 
TACS  is  very  flexible,  to  allow  tailoring  to  a  specific  geographical 
threat,  but  for  this  thesis,  can  be  modeled  as  in  Figure  2: 


Threat 
Prosecut ion 


(More  MCE  units  at  lower  levels  as  required) 


Figure  2.  Typical  TACS  organization 


The  Tactical  Air  Control  Center  (TACC)  is  assigned  overall  control  of 
the  system  and,  under  present  TACS  organization,  has  a  lateral  information 
flow  relationship  with  the  threat  prosecution  assets.  The  Modular  Control 
Elements  (MCE)  are  the  distributed  decision-makers  of  the  system,  and  can 
assume  various  levels  in  the  TACS  hierarchy,  as  depicted  in  Figure  2.  The 
focus  of  this  thesis  will  be  on  the  direct  relationship  between  TACC  and 
one  subordinate  MCE.  The  results  of  this  treatment  will  be  equally  appli¬ 
cable  to  all  subordinate  MCEs  on  an  individual  basis. 


3 .  THE  MODEL 


We  mentioned  in  Chapter  1  that  extensive  work  has  already  been  done  in 

■a 

the  area  of  modeling  the  components  of  C  organizations  and  associated 
information  structures.  Therefore,  no  attempt  is  made  in  this  thesis  to 
provide  a  realistic,  detailed  model.  Instead,  a  very  simplistic  model  is 
developed,  so  that  the  concepts  of  incentive  control  may  be  clearly  pre¬ 
sented  and  easily  understood. 

3.1  General  Structure  of  the  Follower  (MCE) 

The  first  task  in  developing  a  suitable  model  for  the  MCE  (and  ultima¬ 
tely,  its  objective  function)  is  to  distinguish  between  two  subsystems  at 
work  within  MCE: 

(1)  The  command  subsystem,  which  is  MCE  commander  or  the  command 
structure. 

(2)  The  execution  subsystem,  which  carries  out  the  MCE  function,  at 
the  direction  of  the  commander. 

The  command  subsystem,  for  the  purposes  of  this  thesis,  simply  has  the  job 
of  maximizing  the  MCE  objective  function,  J,  given  n  and  determining  a 
value  of  0.  Using  terminology  developed  in  [5],  the  parameter  n  will  be 
the  number  of  "customers”  which  MCE  must  process.  The  output  of  the  MCE 
commander  is  r,  which  we  call  the  "customer  processing  rate."  This  dic¬ 
tated  rate  then  becomes  the  input  to  the  execution  subsystem,  which  is  also 
influenced  by  n  and  0,  but  simply  reacts  to  the  input  r.  Finally,  the  out¬ 
put  of  the  execution  subsystem  is  the  objective  function,  J,  evaluated  at  r, 
It  is  for  TACC,  then,  to  sample  the  output,  r,  of  the  command  sub¬ 
system,  and  observe  J,  to  determine  if  MCE  needs  control.  TACC,  in 


general,  does  not  have  access  to  3,  which  is  MCE  commander's  perception  of 
the  internal  parameter,  but  must  generate  an  estimate,  3'.  Finally,  if 
TACC  does  decide  to  intervene,  the  resulting  n  which  he  computes  is 
announced  to  the  command  subsystem  but  enforced,  if  necessary,  at  the  exe¬ 
cution  sybsystem.  Figure  3  shows  a  block  diagram  of  the  complete  system. 


Figure  3.  Block  diagram  of  TACC-MCE  interaction 
3.2  Modeling  of  MCE 

The  next  step  in  the  process,  as  it  applies  to  TACS,  is  to  model  the 
MCE  execution  subsystem.  One  of  the  more  documented  methods  of  modeling 
human  system  elements  is  the  transfer  function  approach  [6].  It  is  par¬ 
ticularly  suited  to  our  problem  since  r,  the  MCE  independent  rate  variable, 
is  easily  viewed  as  a  frequency  domain  variable  (operations/time  unit). 

The  main  problem  encountered  in  the  transfer  function  approach  to 
human  modeling  is  that  it  is  almost  universally  accepted  that  human  com¬ 
ponents  are  non-linear,  while  the  transfer  function  is  a  linear  analysis 
tool.  The  transition  from  a  non-linear  model  to  a  linear  one  is  provided 


by  [7],  where  linear  and  non-linear  components  of  the  human  system,  in  the 
frequency  domain,  are  combined  as 

J(r)  =•  H(r )  •  X(r)  +  Jr  (4) 

where  J(r)  is  the  output,  H(r)  is  a  linear  transfer  function,  X(r)  is  the 
input,  and  is  a  “remnant"  term,  containing  sytera  non-linear  and  random 
elements.  It  is  incumbent  upon  the  TACS  engineer,  therefore,  to  justify 

assumptions  which  make  negligibly  small.  These  assumptions  must  deal 
with  environmental  factors,  human  inconsistencies  and  just  simply  non¬ 
linear  subtasks.  The  first  two  of  these  factors  can  be  very  subjective, 
while  the  last  can  often  be  handled  by  linear  approximation  techniques. 

The  easiest  way  to  guarantee  linearity  is  to  model  the  subsystem  using 
only  four  basic  sub-elements:  integration,  differentiation,  algebraic 
addition  and  multiplication  by  a  constant  [8].  In  [6],  there  are  several 
possible  combinations  of  these  building  blocks  discussed  for  a  variety  of 
human  operations  In  tracking.  One  of  the  simplest  versions  is  the  classi¬ 
cal,  first-order,  low-pass  filter,  as  shown  in  Figure  4.  Note  that  the 
forward  gain  parameter,  X,  Is  actually  composed  of  two  cascaded  gains,  one 
which  can  be  controlled  externally  and  one  which  cannot.  Obviously,  an 
Increase  in  either  n  or  y  will  result  in  a  decrease  in  forward  gain  of  the 
system.  For  simplicity,  the  feedback  gain  has  been  made  unity,  although 
this  restriction  is  removed  in  the  two-parameter  discussion  in  Appendix  2. 


Figure  4.  Simple  model  of  the  MCE  execution  subsystem 


We  now  go  to  the  frequency  domain,  where  the  model  Is  easier  to  use 
and  understand.  Using  the  variables  as  Indicated  In  Equation  (4)  and 
Figure  4,  we  have: 


X(r) 


*J(r) 


H(r)  -  1/ (jynr  +  1) 


(5) 


Disregarding  the  phase  information,  we  keep  only  |H(r)|.  The  sketch  of 
Figure  5(a)  shows  the  low-pass  characteristics  of  |H(r)|  and  yields  two 
important  (and  logical)  Insights  Into  this  model: 


(1)  Zt  recognizes  that  at  higher  frequencies,  performance,  on  an 
operation  by  operation  basis,  declines. 

(2)  It  Incorporates  both  internal  and  external  parameters  which  can 
greatly  alter  the  transfer  characteristic. 


CAT] 


The  transfer  function  of  Figure  5(a)  is  still  difficult  to  work  with, 
however,  so  it  will  be  approximated  by  a  different  function.  Mote  that  the 


u;  |Hir)|  -  1  /  /  (Ynr)‘  +  1 


(b)  G(r)  -  2  -  ep 


Figure  5.  Transfer  characteristics  of  the  MCE  model 

convex  portion  of  |H(r) | ,  in  Figure  5(a),  is  similar  to  an  Inverted  expo¬ 
nential.  In  fact,  if  |H(r)|  is  truncated  at  rc,  it  can  be  approximated  by 

G(r)  -  2  -  e0nr,  (6) 

where  S  *  ky,  for  some  constant  k.  In  this  form,  we  have  the  familiar 
Internal  parameter,  0.  Figure  5(b)  demonstrates  how  closely  the  truncated 
portion  of  |H(r)|  may  be  approximated  by  G(r). 

The  final  step  in  the  modeling  process  is  to  form  J(r),  the  objective 
function  of  MCE  and  the  output  of  the  MCE  execution  subsystem.  Recall  from 
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Section  2.1  that  the  input  to  the  execution  subsystem  was  simply  r,  the 
operating  rate.  Starting  with  the  general  form  of  J(r),  we  have: 

J(r)  =*  |H(r)|  •  |x(  r )  | 

J(r)  =  G(r)  •  |X(r) | 

-  r  (2  -  e0™")  (7) 


Equation  (7)  is  the  final  form  of  the  model,  which  will  be  used  for  the 
remainder  of  the  thesis  (except  Appendix  2).  It  is  convex  for  all  positive 
values  of  r,  n,  and  0,  and,  if  n  is  treated  as  a  constant  (no  control).  Is 
maximized  as  follows: 


J(r) 


max 


J(r*) 


r*  -  {r:  dJ(r)/dr  -  0} 


dJ/dr  -  2  -  (0nr*  +  1)  e^nr*^  *  0 
r-r* 


r*  -  .375/fJn 

J(r*)  -  (,375/0n)  (2  -  e(*375))  -  .2O4/0n 


(8) 

(9) 


Figure  6  compares  the  objective  function,  J(r),  for  different  values  of  n, 
with  0  fixed,  and  different  values  of  0,  with  n  fixed.  It  also  depicts  the 


location  of  r*  for  each  curve 


Figure  6.  The  ideal  MCE  objective  function 


In  summary,  we  have  developed  a  model  for  the  execution  subsystem  of 
MCE,  which  is  the  frequency  domain  objective  function  J(r).  The  input,  r, 
is  a  rate  value  dictated  by  the  command  subsystem,  which  must  maximize 
J(r).  The  model  contains  two  parameters:  8,  an  internal  parameter,  which 
is  characteristic  of  the  execution  subsystem;  and  n,  the  number  of  custom¬ 
ers,  which  can  be  affected  externally.  Ve  now  turn  our  attention  to  the 
manner  in  which  this  model  may  be  controlled. 


4.  THE  PIECE-WISE  LINEAR  APPROACH 


A  significant  problem  In  the  incentive  control  approach  is  finding  the 
strategy  which  will  bring  about  the  desired  response  in  the  follower.  In 
[3]  and  [4]  it  is  seen  that  if  the  strategy  is  made  a  function  of  the  dif¬ 
ference  between  the  actual  and  desired  responses,  the  follower,  when  he 
substitutes  the  incentive  into  his  own  objective  function,  will  be  forced 
to  maximize  at  or  near  the  desired  point. 

In  Section  2.1,  the  variable  n  became  the  leader-follower  control 

A 

variable,  as  a  function  n.  This  is  the  same  n  as  in  Equation  (7)  of  the 
MCE  model,  and  it  will  be  the  control  for  MCE,  with  TACC,  the  leader, 

A 

determining  n.  It  is  tempting,  at  first,  to  use  a  linear  construction  of 
the  control: 


n  *  nQ  +  k  (r  -  r*),  nQ  >_  0,  k  _>  0  (10) 

where  n  Is  a  constant  term,  representing  the  number  of  customers  MCE  will 
o 

have  to  handle  no  matter  what  (TACC  may  have  no  control  over  nQ);  r  is  the 
Independent  rate  variable;  r*  is  the  solution  to  Equation  (8),  given  nQ  and 
S';  and  k  is  a  gain  parameter  which  is  calculated  by  TACC. 

The  problem  with  Equation  (10)  is  that  if  (r  -  r*)  <  0,  then  n  <  nQ. 

It  is  seen  from  Figure  6(a)  that  smaller  and  smaller  values  of  n  lead  to 
larger  and  larger  maximum  values  of  J(r)  and  allow  higher  r*.  In  fact,  if 
(r  -  r*)  ■  -nQ,  then  n  -  0  and  MCE  has  no  customers  to  process  and  achieves 
its  highest  maximum!  It  was  stipulated,  however,  that  nQ  was  the  minimum 
number  of  customers  that  must  be  serviced,  so  we  are  led  to  propose  a 
slightly  different  n. 


4.1  The  Control  Structure 


Consider  the  following  control  structure: 


n  =  n  +  k  r  -  r*  . 

o  1  ' 


With  this  function,  which  shall  be  called  the  piece-wise  linear  incentive 
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strategy,  n  >_  nQ  for  any  r  and  k  _>  0.  In  view  of  Figure  6(a),  MCE  can  do 
no  better  than  processing  nQ  customers  at  r*: 

n|  .an+k|r*-r*|=n  (IS 

r»r*  o  1  1  o 


Before  the  effect  of  this  control  is  investigated,  it  must  be  noted 
that  the  piece-wise  linear  approach  does  not  lead  to  a  continuously  dif¬ 
ferentiable  objective  function  in  r,  but  contains  a  discontinuous  slope  at 
r  -  r*  (see  Figure  7(c)).  Therefore,  Equation  (3)  cannot  be  used  to  find 
r*.  Instead,  the  slope  of  the  objective  function  must  be  calculated  on 
each  side  of  the  "peak"  at  r*,  to  ensure  that  the  largest  value  of  J(r) 
occurs  at  that  peak. 

Suppose  TACC  announces  a  strategy  of  the  form  in  Equation  (11).  The 
MCE  commander,  in  order  to  effectively  maximize  his  objective  function, 
must  substitute  n  for  n  in  Equation  (7): 


J(r)  -  r  (2  -  e0nr) 


-  r  (2  -  e0[n°  +  klr-r*U  r) 


Now  equation  (13)  must  be  analyzed  for  two  intervals  of  the  r-axis:  r  ^  r* 
and  0  <  r  <  r*. 


Statement. :  Suppose  nQis  known  and  0  '  is  estimated  by  TACC,  so  that  r* 


can  be  calculated  by  Equation  (8).  Suppose  also  that  the  minimum 


expected  value  of  8,  denoted  0 ^ ,  is  known.  Then  a  sufficient  con¬ 


dition  for  J(r)  to  have  a  negative  slope  on  the  interval  [r*,°»]  is  to 
choose  a  k  satisfying: 

(0 ,n_r*) 

(14) 


k  > 


2  -  (0 .nr*  +  1)  e  1  ° 
l  o 


2  (8 .nr*) 

0Lr*^  e  1  ° 


Proof :  The  objective  function,  for  this  interval,  can  be  written: 


J(r)  -  r  (2  -  e9ln°  *  k(r'r*,lr) 


(15) 


If  J(r)  has  a  negative  slope  for  all  r  in  [r*,«]  and  8  in  a  given  range, 
8  >  8j,  then  its  maximum  value  will  occur  at  r*.  This  can  be  stated 
mathematically: 


dJ/dr  <  0 

A 

2  -  [0  (n  +  dn/dr)  r  +  l]e^nr^  <  0 

A 

dn/dr  -  k,  based  on  Equation  (15), 


so  2  -  (0nr  +  1)  e(Snr)  -  k0r2  e(3nr)  <  0 


(16) 


The  worst  case  for  ensuring  Inequality  (16)  is  if  0  and  r  are  at  their 
minimum  values.  Therefore,  let  0  *  0^,  r  ■  r*,  and  thus  n  *  nQ. 
Inequality  (16)  becomes: 


vv. 


Therefore, 


(0 ,n  r*) 

2  -(8  n  r*  +  1)  e  1  0 
_ 1  o _ 

2  (0 ,n  r*) 

8,r*  e  1  o 


Therefore,  when  TACC  announces  this  strategy,  he  will  have  to  select  k 
based  on  Inequality  (14).  However,  there  is  another  side  of  the  function 
which  must  be  considered. 

4.1.2  Case  II:  0  <  r  <  r* 

This  case  is  a  little  more  difficult  to  analyze.  Normally,  it  would 
be  sufficient  to  guarantee  a  positive  slope  at  r*,  but  if  k  is  made  large 
enough,  a  cubic-like  oscillation  (maximum  followed  by  a  minimum)  is  induced 
on  [0,r*].  Therefore,  k  must  be  large  enough  to  provide  a  positive  slope 
at  r*,  but  not  so  large  as  to  induce  an  oscillation  between  0  and  r*. 

Statement:  Suppose  nQ  is  known  and  0’  estimated  by  TACC,  so  r*  can  be 
calculated.  Suppose  that  the  maximum  expected  value  of  0,  denoted  0£, 
is  also  known.  Then  a  sufficient  condition  for  J(r)  to  have  a  positive 
slope  on  the  interval  [0,r*]  is  to  choose  a  k  satisfying: 

(0~n  r*) 

2  -  (0„n  r*  +  1)  e  2  ° 


2  (0,n  r*) 

02r*^  e  2  ° 

provided  there  is  no  real  solution,  in  the  interval  [0,r*],  to  the 
equation  dJ/dr  *  0. 

Proof :  For  this  interval,  the  objective  function  becomes: 

0 [n  +  k(r*-r) ]r 

J(r)  »  r(2  -  e  °  )  ( 


As  mentioned  earlier,  the  first  task  is  to  guarantee  a  positive  slope  at 
r  -  r*,  with  0  in  a  given  range  (0  <  0„): 
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dJ/r  >  0  and  dn/dr  =  -k 


2  -  (gnr  +  1)  e(Snr)  +  kgr2  e(f3nr)  >  0. 


Finally,  making  the  substitutions  r  =  r*,  g  =  and  n  =  n^: 

(0„n  r*)  2  (8,n  r*)  >  0 

2  -  (g„n  r*  +  1)  e  2  °  +  kg„r*^  e  2  0 


Therefore, 


k  >  - 


(g„n  r*) 

2  -  (0,n  r*  +  1)  e  2  ° 
_ l  o _ 

2  (8,n  r*) 

&2r*  e  2  o 


In  summary,  if  TACC  decides  to  use  the  piece-wise  linear  approach,  he 
must  calculate  a  minimum  value  of  k  for  both  sides  of  the  slope  discon¬ 
tinuity  at  r*.  He  then  must  select  a  k  which  satisfies  both  cases  but  is 
not  so  large  that  a  maximum-minimum  pair  forms  on  [0,r*].  The  following 
example  will  help  greatly  in  understanding  the  preceding  discussion. 


4.2  An  Example 

Suppose  that  TACC  estimates  that  MCE  should  have  an  internal  parameter 
8'  *  .01,  and  it  is  given  that  nQ  *  10.  Figure  7(a)  depicts  the  ideal  MCE 
objective  function  (without  control).  The  r*  corresponding  to  the  maximum 
is  found  by  Equation  (8): 


r*  -  .375/8  'n  -  . 375/ ( . 0 1 ) ( 1 0 )  -  3.75 
o 


It  is  assumed  that  TACC's  estimate  of  6 '  is  correct,  but  that  MCE  commander 
perceives  a  different  internal  parameter,  8  ■  .008.  Therefore,  MCE  will 
use  a  different  objective  function  to  maximize,  which  is  depicted  in  Figure 
7(b).  Because  of  this,  he  will  actually  drive  the  MCE  at  some  r  *  r*, 


2 


v.  <r 
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where 


r  =  .375/ ( .008) (10)  =  4.69 


When  TACC  observes  MCE  operating  at  this  unexpected  rate,  he  may  wish  to 
intervene,  using  the  piece-wise  linear  approach. 

First,  TACC  must  estimate  the  range  in  which  6  falls.  Suppose  the 
estimated  range  is  .006  <  6  <  .012. 


Case  I:  r  >  3.75 


(8,n  r*) 

2  -  (fin  r*  +  1)  e  1  0 
1  o 


2  (8 ,n  r*) 

tf*  e  1  o 

and  if  8^  *  .006,  nQ  =  10,  and  r*  =  3.75,  then  k  >  4.41 


Case  II:  0  <  r  <  3.75 


k  >  - 


(8,n  r*) 

2  -  (g  n  r*  +  1)  e  2  0 
_ z  o _ 

2  (89n  r*) 

82r*  e  2  0 


and  82  =  -012,  so  k  >  1.04. 

Therefore,  in  order  to  satisfy  both  Inequalities  (20)  and  (21),  TACC 
decides  to  pick  k  =  5.  The  incentive  strategy  that  he  announces  is: 


n  =  10+  5  r  -  3.75  . 


Once  MCE  substitutes  Equation  (22)  into  his  objective  function,  his  new 
curve  will  appear  as  in  Figure  7(c).  Note  that  the  function  must  be  maxi¬ 
mized  at  r  *  r*  ■  3.75  and  that  no  cubic  oscillations  occur  on  [0,r*]. 
Figure  7(d),  however,  demonstrates  how  a  large  enough  gain  (k  =  17)  can 
induce  the  maximum-minimum  pair. 


■  .v.v.'A 


1.0  2.0  3.0  4.0  5.0  r 


(c)  With  control:  k  ■  5 


L.U  ^.U  J.U  4 


(d)  With  control:  k  ■  1 


Figure  7.  Objective  function  with  piece-wise  linear  control 


It  is  important  to  remember  that  although  8  =  .008  is  used  to  generate 

curves  (b),  (c),  and  (d)  in  Figure  7,  any  8  in  the  interval  [.006,. 012] 

would  have  yielded  the  same  r*,  although  the  shapes  of  the  curves  and  the 

resulting  J  values  would  vary.  Because  of  this,  r*  is  considered  insen— 
max  ’  - 

sitive  to  changes  in  3,  for  the  range  .006  <  3  <  .012. 


5.  THE  QUADRATIC  APPROACH 


5.1  The  Control  Structure 

There  is  an  alternative  approach  to  controlling  the  MCE,  and  that  is 
to  use  a  control  variable  of  the  form: 

n  =  n  +  k.  (r  -  r*)2  +  k_  (r  -  r*)  (23) 

oi  i 

In  this  case,  nQ  and  r*  are  the  same  as  in  the  piece-wise  linear  approach, 
but  now  two  gain  constants,  k^  and  k2,  have  been  introduced.  One  of  the 

a 

advantages  of  this  structure  is  that  n  is  a  continuous  function,  so  MCE’s 
objective  function  will  also  be  continuous  when  this  control  is  substituted 
for  n.  Therefore,  J(r)  does  not  have  to  be  evaluated  for  two  different 
cases,  as  in  the  former  approach,  and  r*  can  be  found  using  Equation  (3): 

r*  -  {r:  dJ/dr  »  0}  (24) 

dJ/dr  ■  2  -  |8r[(dn/dr)r  +  n]  +  l}  e^nr)  ■  0  (25) 

where  dn/dr  *  2k^(r  -  r*)  +  k2>  Therefore,  by  substitution: 

dJ/dr  *  2  -  {f)r[(2kj(r-r*)  +  k2)r  +  (nQ  + 

i  r  .  y.  r  *01  .  ii  (8nr)  -  0  (26) 

k^r-r*)  +  k2(r-r*)J]  +  1}  e 

The  solution  of  Equation  (26)  is  not  a  trivial  matter,  and  we  found  it 

convenient  (and  very  accurate)  to  use  the  following  procedure: 

dJ/dr  -  2  -  (f2(r)  +  f^r)  +  1]  e(Snr)  -  0 

where  f^(r)  and  f2(r)  are  determined  from  Equation  (26).  Thus 

r  -  { ln{ 2 )  -  ln[f2(r)  +  fj(r)  +  1)}/Sn 


Using  a  first  order  Taylor  series  approximation  of  the  second  logarithmic 
term  about  r*,  the  problem  eventually  reduces  to  a  cubic  equation,  which  is 
solved  by  computer  using  the  cubic  formula. 

Calculation  of  the  constants  kj  and  k2  by  TACC  is  relatively  easy, 
however,  if  we  remember  that  TACC's  goal  is  to  have  Equation  (26)  satisfied 
at  r=r*.  By  substituting  r*  into  Equation  (26),  it  is  seen  that  most  terms 
drop  out,  and  the  remaining  equation  is: 


(Bn  r*) 

2  -  (Br*[k2r*  +  nQ]  +  1)  e  °  -  0 


(Bn  r*) 

2  -  (BnQr*  +  1)  e  ° 

2  (Bn  r*) 

Br*  e  ° 


It  is  interesting  to  note  that  Equation  (27)  depends  only  on  k^.  This  is 
significant  in  that  if  TACC  knew  exactly  which  B  MCE  was  going  to  assume  in 
its  maximization,  he  could  simply  calculate  k2  by  Equation  (27)  and  use  any 
kj,  with  the  resulting  objective  function  maximized  at  r*. 

In  reality,  however,  B  is  not  known  by  TACC,  and  it  is  desired  to 
minimize  the  sensitivity  of  r*  to  fluctuations  in  B.  Again,  it  is  assumed 
that  B  can  be  confined  to  a  range,  such  as  0^  <_  B  .<  Additionally,  an 

average  B  for  the  range,  Bc»  must  be  selected  by  TACC.  Equipped  with  this 
data,  TACC  can  methodically  select  k^  and  k2  so  that  if  MCE  uses  Bc,  he 
will  maximize  exactly  at  r*,  and  if  he  uses  some  other  B  in  [0p02],  his 
maximization  will  remain  within  a  given  tolerance  of  r*.  One  additional 
advantage  of  this  policy  is  that  B’»  the  value  which  TACC  thinks  actually 
characterizes  MCE  and  upon  which  r*  is  based,  does  not  have  to  be  in  the 
Interval  [Bp  Bj]*  ^  0'  does  fall  in  [0p02],  however,  it  is  convenient 
to  let  8C  *  B'  for  analysis  and  calculation  purposes. 
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The  procedure  for  finding  and  k£  is  as  follows: 

(1)  Define  the  anticipated  range  of  3,  by  selecting  3j,  3,  and  3^. 

(2)  Compute  r*  for  the  ideal  objective  function  and  specify  the 

allowable  tolerance  in  r,  the  operating  rate  of  MCE,  about  r*. 

This  is  also  expressed  as  an  interval,  [r^,^]'  containing  r*. 

(3)  Using  n^,  r*  and  0c>  calculate  k?  using  Equation  (27). 

(4)  Using  equation  (26)  and  the  above  computed  value  of  k^ ,  find  kj 
at  each  boundary  of  (rj.^J: 

a.  At  r^ ,  use  3 ^ 

b.  At  r- ,  use  0  ^ 

The  final  k^  that  TACC  selects  must  be  larger  than  the  largest 
kj  computed  for  both  of  the  above  cases. 

There  is  one  more  consideration  in  selecting  k^:  if  k^  is  too  large, 
the  induced  MCE  objective  function  will  have  a  very  sharp  "peak"  at  its 
maximum.  This  could  mean  that  even  if  MCE  is  maximizing  within  [rj,^] 
tolerance,  the  maximization  point,  r,  may  still  be  far  enough  away  from  r* 
to  cause  a  substantial  increase  in  n.  Recall  that  n  is  based  on  deviations 
in  r  from  the  ideal  r*.  If  k^  is  large  enough,  even  a  small  deviation  in  r 
will  cause  n  to  be  enforced  (i.e.,  extra  customers  are  given  to  MCE  for  not 
"cooperating").  This  is  inconsistent,  however,  with  the  fact  that  MCE  is 
within  tolerance. 

The  solution  to  this  dilemma  lies  in  analyzing  the  "discrete"  nature 
of  n.  Since  TACC  cannot  give  fractions  of  customers  to  MCE,  only  Integral 
values  of  n  are  permitted.  Therefore,  TACC  must  wait  until  n  increases 
sufficiently  over  nQ  to  actually  increment  with  an  additional  customer 
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(enforce).  For  this  thesis,  we  assume  that  if  n  *  +  .5,  the  control 

will  be  enforced.  This  discussion  leads  us  to  a  fifth  step: 

(5)  Compute  kj  ,  such  that  for  r  in  [r1,r2J,  n  <  n  +  .5.  This 
max  ° 

last  step  may  be  accomplished  precisely  by  evaluating  Equation 

(23)  at  r1  and  or  by  using  a  graph  of  n  vs.  r,  as  in  Figure 

8(c). 


5.2  An  Example 

Again,  suppose  nQ  -  10  customers  and  TACC  uses  0'  *  .01  as  the  actual 
parameter  of  MCE.  It  is  anticipated,  however,  that  MCE  will  choose  a  0  in 
the  interval  .007  £  8  £  .015.  For  simplicity,  TACC  uses  0c  -  .01.  Suppose 
further  that  the  maximum  allowable  tolerance  in  r  is  3.6  r  3.9  (recall 
that  r*  for  these  data  is  3.75).  Proceeding  from  the  third  step: 


(3) 


(8  n  r*) 

2  -  (8  nr*  +  1)  e  c  ° 
c  o 

2  (8  nr*) 

0  r*  e^  c  ° 
c 


2  -  [(.01)(10)(3.75)  +  1]  (1.455) 
(.01)(3.75)2(1.455) 


(4)  At  Tj  ■  3.6,  02  ■  .015.  Since  kj  *  0,  using  Equation  (26): 

kj  -  7.1 

At  ■  3.9,  m  *007.  With  kj  •  0  and  Equation  (26): 

kj  -  7.5 

Therefore,  must  be  greater  than  7.5. 


0,  Equation  (25)  becomes: 


(5)  Since  k2  * 

n  -  10  +  kj ( r  -  r* )“ 

Evaluating  at  r^  *  3.6: 

n  =  10  +  kj(3.6  -  3. 7 5) 2  <  10.5 

kj  <  22.2 

Because  k2  ■  0,  Equation  (23)  is  symmetric  about  r*  and  the  same 

result  is  obtained  for  r2  ■  3.9.  Therefore  k^  ■  22.2. 

max 

Thus,  TACC  must  select  k2  ■  0  and  k^  in  the  range  [7.5,  22.2}.  Suppose  he 
chooses  kj  -  10.  The  control  that  is  announced  by  TACC  to  MCE  is: 

n  «  10  +  10(r  -  3.75)2  (28 

Finally,  for  this  example,  assume  that  MCE  is  actually  using  8  «  .008  in 
its  operation.  After  substitution  of  Equation  (28)  in  the  MCE  objective 
function,  the  MCE  must  maximize: 

j,  _  r(2  _  e.008[10  +  10(r  -  3.75)2]r j  (29 

Figure  8  is  an  illustration  of  this  example.  In  Figure  8(a),  the  MCE 
objective  is  depicted  for  8  *  .008  without  control.  Note  that  the  maximum 
occurs  at  r  ■  4.69.  In  Figure  8(b),  the  effect  of  the  control  is  seen  as 
Equation  (29)  is  displayed.  MCE  maximizes  at  r  *  3.82,  which  is  well 
within  tolerance.  Figure  8(c)  is  a  plot  of  n  vs.  r,  about  r*  ■  3.75,  for 
different  values  of  kj  (k2  ■  0).  The  superimposed  rectangle  is  the  region 

A 

in  which  n  is  essentially  nQ,  because  of  the  discrete  nature  of  the 
control.  Figure  8(c)  may  serve  as  a  "look-up  table"  for  determining  the 


maximum  value  of  kj  (step  5  in  the  procedure).  The  purpose  of  Figure  8(d) 
is  to  demonstrate  the  "centering"  properties  of  k2«  The  data  are  the  same 
as  for  Figure  8(b),  except  the  objective  function  is  sketched  for  k2  “  5 
and  k2  -  -5.  If  TACC  knew  that  8  -  .008  in  advance,  the  actual  centering 
value  of  k2  could,  of  course,  be  calculated  from  Equation  (27). 

2  -  [  ( .008)  ( 10)  (3.7  5)  +  1]  (1.35)  _  u6n 
2  (. 008)(3. 75)2(1. 35) 


1.0  2.0  3.0  4.0  5.0  r 


l.CT  2.0  3.0  4.0  r 


(a)  No  controls  0  ■  .008 


(b)  With  control:  k,  -  10,  k. 


Figure  8.  Objective  function  with  quadratic  control 


6.  CONCLUSION 


In  this  thesis,  we  have  developed  the  mathematical  framework  for 
employing  incentive  control  in  the  Tactical  Air  Control  System.  We  modeled 
MCE  as  being  composed  of  a  command  subsystem,  whose  function  is  to  maximize 
its  objective  function,  and  an  execution  subsystem,  which  behaves  like  a 
transfer  function,  whose  characteristics  are  governed  by  an  internal  and  an 
external  parameter. 

The  goal  of  the  Incentive  control  strategies,  introduced  in  Chapters  4 
and  5,  is  to  ensure  that  MCE  maximizes  its  objective  function  at  or  near  a 
desired  operating  point,  r*,  regardless  of  the  value  of  the  Internal  para¬ 
meter  MCE  uses  in  its  objective  function.  The  piece-wise  linear  control 
strategy  makes  r*  insensitive  to  changes  in  0  over  a  predetermined  range 
of  0,  but  is  computationally  more  difficult  to  analyze.  The  quadratic 
approach  is  somewhat  easier  to  analyze  and  has  the  advantage  that  0 '  does 
not  have  to  be  in  [0^,02].  It  is  not,  however,  totally  insensitive  to 
changes  in  0  but  can  be  made  minimally  sensitive  for  given  tolerances  in 
MCE  response.  Both  controls  have  the  distinct  advantage  that  for  any  model 
selected,  if  the  objective  function  is  a  monotonically  decreasing  function 
of  the  control  variable  n  (reference  Figure  6),  the  induced  objective  func¬ 
tion  will  be  maximized  only  at  the  bias  value,  nQ  (which  corresponds  to 
r*). 

3 

Another  important  application  of  incentive  mechanisms  in  the  C 
setting  is  in  Command,  Control,  and  Communications  Countermeasures  (C  CM). 
In  this  context,  the  participants  are  viewed  as  non-cooperative.  The  use 
of  incentive  strategies  by  one  superpower  to  induce  "cooperation"  by 
another  superpower  in  strategic  arms  control  is  discussed  in  [4].  These 
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same  ideas  can  be  extended  to  the  tactical  theater  to  force  an  enemy 
C  system  to  behave  in  a  certain  manner.  The  steps  taken  are  exactly  the 
same  as  outlined  in  this  thesis: 

(1)  Model  the  adversary's  operation  and  determine  his  objective 
function. 

(2)  Determine  all  unknown  internal  parameters  of  the  objective  func¬ 
tion. 

(3)  Determine  what  parameter(s)  in  the  objective  function  can  be 
affected  by  an  incentive  strategy. 

(4)  Design  an  appropriate  strategy. 

(5)  Announce  the  strategy  to  the  adversary. 


APPENDIX  1 


SENSITIVITY  FUNCTION  ANALYSIS 


Throughout  this  thesis,  our  discussion  has  been  based  on  the  fact  that 

MCE  might  choose  some  value  of  6  other  than  8'  in  its  objective  function. 

At  the  same  time,  our  assumption  has  been  that  8'  (the  value  that  TACC 

thinks  MCE  should  have)  is  actually  the  more  appropriate  value.  We  desire, 

then,  to  have  some  measure  of  how  close  J  ,  the  actual  maximum  value  of 
’  max  - 

the  MCE  output,  is  to  J’  ,  the  ideal  calculated  by  TACC. 

max 

In  [4],  sensitivity  functions  are  introduced  to  provide  such  a 
measure.  Using  the  derivative  terms  of  the  Taylor  series  expansion  of  J 
about  8',  when  r  =  r*,  we  denote  the  first  order  sensitivity  function: 


I,  -  dJ(n,r,8)/d8 
1  8-8* 
r=r* 
n=n 

°  (30) 

=  [ (dJ/dn) (dn/ dr)  +  dJ/dr]  (dr/d8) 

8=8* 

r»r* 

n*n 

o 


We  first  analyze  each  term  of  Equation  (30)  in  view  of  the  piece-wise 
linear  approach.  dJ/dn  will  always  be  negative,  due  to  the  monotonic  rela¬ 
tionship  between  J  and  n  in  the  original  model,  dn/dr  is  discontinuous  at 
r  *  r*,  since  it  is  equal  to  k,  for  r  >  r*,  and  -k,  for  r  <  r*.  dr/d8  is 
zero,  however,  because  in  the  piece-wise  approach,  we  have  total  insen¬ 
sitivity  of  r  to  changes  in  8  in  [0  ^ ,0  £  J  *  which  contains  8'.  Due  to  this 
fact,  the  entire  function  is  zero. 

For  the  quadratic  approch,  we  can  never  make  dr/d8  equal  to  zero  but 
can  only  minimize  it.  dJ/dr  is  still  zero,  of  course,  and  dJ/dn  is  still 
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always  negative.  The  easiest  way  to  make  1^  =  0,  then,  is  to  ensure  dn/dr 
=  0.  Recalling  the  discrete  nature  of  n,  we  can  say  Lhat  for  a  s;nalL  range 
of  r  about  r*,  the  change  in  n  with  respect  to  r  is  zero.  Furthermore, 
since  our  selection  of  was  based  on  a  range  of  3 ,  if  3 '  is  included  in 

^1*^2^*  t*ien  ^n/dr  =  0  over  that  range  (in  a  "discrete"  sense).  If  3  '  is 
not  included  in  we  cannot  use  this  function  as  a  measure  of  sen¬ 

sitivity  about  3'. 

If  1^  =  0,  the  dominant  term  of  the  expansion  becomes  the  second  order 
sensitivity  function,  which  we  define  as: 

I,  =  d2J(n,r,8)/d82  I 


2  2 

=>r0(nJ  +  n  J  +nJ  +J) 

8  r  nn  r  nr  r  rn  rr 


+  (J  n  +  J  )r00 
n  r  r  88 


Again  analyzing  the  piece-wise  linear  approach  first,  we  find  that  since 


dr/d8  *  0  at  8  “  8',  the  entire  first  term  is  zero.  Total  insensitivity  of 
r  to  8  in  C0 1  *6 2 3  also  guarantees  that  r^  »  0,  thus  I2  *  0  .for  the  piece- 
wise  approach.  In  fact,  all  higher  order  derivatives  of  J  will  also  vanish 
in  a  neighborhood  of  8'»  so  for  this  approach  the  maximum  value  of  the 
objective  function  is  totally  Insensitive  to  changes  in  8  in  [8^,82]* 

For  the  quadratic  approach,  since  dJ/dr  =■  0  and  dn/dr  ■  0,  Equation 
(31)  reduces  to  the  following  expression: 


I2  -  (dr/d8)2(d2J/dr2) 


(32) 


We  know  that  dr/d0  cannot  be  zero  but  is  minimized  for  kj  very  large. 
Therefore,  the  first  product  term  can  only  be  a  small  positive  number. 
Also,  the  objective  function  is  convex  in  r  for  all  3  and  n,  so  the  second 
product  terra  is  always  a  negative  valued  function.  Therefore,  to  minimize 


I 2»  we  should  make  kj  as  large  as  possible,  but  not  larger  than  kj 
calculated  in  Chapter  5  (Equation  27). 


max 
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APPENDIX  2 

A  TWO-PARAMETER  MODEL 

In  Section  3.2,  a  model  was  developed  for  MCE,  which  included  one 
internal  parameter,  0.  Now  we  investigate  the  effect  of  incorporating  a 
second  parameter,  a,  into  the  objective  function.  Consider  a  transfer 
function  of  the  form: 

G(r )  =  a  -  e(enr)  (33) 

The  difference  between  Equation  (33)  and  the  one-parameter  case  corresponds 
roughly  to  changing  the  unity  feedback  gain  of  Figure  4  to  a  variable  feed¬ 
back.  The  effect  of  introducing  a,  for  a  _<  2,  is  to  vary  the  overall 
amplitude  of  the  transfer  function.  The  resulting  objective  function  for 
MCE  is: 


J(r )  -  r(a  -  e(f?nr))  (34) 

Examining  Equation  (34)  for  the  piece-wise  linear  approach  poses  no 
special  problems.  We  again  divide  the  r-axis  into  two  regions  and  compute 
dJ/dr  for  two  cases. 

Case  I:  r  2i  r* 

Rewriting  Equation  (18)  for  the  more  general  case  involving  a: 

a  -  (0nr  +  l)e(enr)  -  k8r2e(5nr)  <  0  (35) 


Now,  in  addition  to  the  worst  case  values  of  8,r,  and  n,  the  largest  value 
of  a  also  must  be  selected,  to  ensure  Inequality  (35).  Thus,  we  use  a  *  2 
The  result  for  this  case,  then,  is  the  same  as  for  the  one-parameter  case: 


The  same  procedure  is  followed  here  as  in  Section  4.1.2,  except  that 

the  minimum  anticipated  value  of  a  is  used  in  calculation  of  k.  Denoting 

this  value  of  a  as  a  : 

m 


k  > 


(S0n  r*) 

a  -  (B0n  r*  +  l)e  2  o 
m _ 2  o _ 

2  (B,n  r*) 

02r*  e  2  ° 


(37) 


The  problem  of  Inducing  the  extra  oscillation,  for  large  k  on  [0,r*],  is 
even  more  acute  for  smaller  values  of  a.  In  fact,  for  a  <  1.5,  there  may 
not  exist  a  k  which  satisfies  Inequality  (37)  but  which  avoids  the  maximum- 
minimum  phenomenom.  Therefore,  for  this  model,  a  should  be  restricted  to 
the  interval  [1.5, 2.0]. 

Analysis  of  the  two-parameter  case  is  likewise  straightforward  for 
the  quadratic  approach.  The  steps  taken  are  identical  to  the  procedure 
outlined  in  Chapter  5,  except  that  an  anticipated  a  range,  ^.ac^.a2’  *s 
Identified  as  well  as  the  f$  range.  Then,  the  ’’centering”  property  of  k2  at 
ac,  &c  is  retained  and  kj  is  selected  for  a  predetermined  r  tolerance. 


evaluated  at  the  boundaries  r. ,  a.,  and  r0,  o~,  £5 ,  • 
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